COSMOLOGY WITH A SHOCK- WAVE 

November 9, 1998 



Joel Smoller 1 Blake Temple^ 



Abstract 

We construct the simplest solution of the Einstein equations that 
incorporates a shock-wave into a standard Friedmann-Robertson- 
Walker metric whose equation of state accounts for the Hubble con- 
stant and the microwave background radiation temperature. This 
produces a new solution of the Einstein equations from which we are 
able to derive estimates for the shock position at present time. We 
show that the distance from the shock-wave to the center of the ex- 
plosion at present time is comparable to the Hubble distance. We are 
motivated by the idea that the expansion of the universe as measured 
by the Hubble constant might be accounted for by an event more sim- 
ilar to a classical explosion than by the well-accepted scenario of the 
Big Bang. 



1 Introduction 

In the standard model for cosmology it is assumed that, on the 
largest scale, the entire universe is expanding at a rate measured 
by the Hubble law, [0, @, |9|, 0, [18], [Hj . Hubble's Law correlates 
recessional velocities of galaxies with red-shifts. However, this 
correlation has only been verified for nearby galaxies, and it is an 
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extrapolation to apply this law to the entire universe. Moreover, 
it follows from the Einstein equations, [f|, |5], |I8j| , that if the uni- 
verse is everywhere expanding, then every spacetime point can 
be traced back to a singularity in the past, a singularity from 
which the entire universe burst in an event referred to as the Big 
Bang. The assumption that the Hubble Law applies to the entire 
universe is referred to as the Cosmological Principle, [|TTf. The 
Cosmolgical Principle is the starting assumption in the modern 
theory of cosmology, and it is precisely this principle that forces 
the singularity into the standard Big Bang interpretation of the 
origin of the universe. In this paper we explore the possibility 
that Hubble's Law actually only measures a localized expansion 
of the universe, and not the expansion of the entire universe. We 
demonstrate the consistency of this possibility by constructing 
the simplest possible solution of the Einstein equations that ac- 
counts for the observed Hubble expansion rate and the correct 
microwave background radiation temperature, such that there 
is a shock-wave present at the leading edge of the expansion. 

Our motivation is the idea that the expansion of the universe, 
as measured by the Hubble constant, might be the result of a 
large scale localized explosion that generated a shock-wave at 
the leading edge, not unlike a classical explosion into a static 
background, except on an enormously large scale. If this were 
true, then it would place our solar system in a special posi- 
tion relative to the center of the explosion, and this would vi- 
olate the so-called Copernican Principle, at least on the scale 
at which the Hubble Law applies. The Copernican Principle 
is the statement that the earth is not in a "special place" in 
the universe. This principle justifies the standard cosmology 
based on the Friedmann- Robertson- Walker (FRW) metric be- 
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cause the FRW metric is the unique metric that is consistent 
with the Einstein equations, and is homogeneous and isotropic 
about every point. The high degree of uniformity of the back- 
ground microwave radiation in all directions, together with the 
directional independence of the redshifting of galaxies, provides 
the strongest support for the Copernican Principle. The idea 
that there is a shock-wave present at the leading edge of that 
portion of the universe where the Hubble constant applies, also 
violates another basic tenet of modern cosmology; namely, that 
we can meaningfully time reverse the continuum model all the 
way back to microseconds after the Big Bang. (For example, 
the theory of inflation applies to the regime 10 -34 to 10 -32 sec- 
onds after the Big Bang in the continuum model, [j8j].) Indeed, 
it follows from the mathematical theory that shock-waves intro- 
duce a fundamental increase of entropy and consequent loss of 
information, [j7], [T7j]. Thus, when a shock-wave is incorporated 
into cosmology, it becomes impossible to reconstruct the details 
of the early explosion from present data, at least at the level of 
the continuum model. 

The simplest shock-wave model for cosmology is one in which 
the "expanding universe", inside the shock-wave, is modeled 
by the standard FRW metric of cosmology, and the spacetime 
on the outside is modeled by a Tolman-Oppenheimer-Volkoff 
(TOV) metric, (the general relativistic version of a static fluid 
sphere), such that the interface in between is an exact, spheri- 
cally symmetric shock-wave solution of the Einstein equations, 
that propagates outward. The assumption that outside the 
shock-wave is a time-independent spherically symmetric solu- 
tion is not unreasonable if one imagines that the spacetime be- 
fore the explosion occured took a long time getting into the pre- 
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explosion configuration; and the assumption that an expand- 
ing FRW metric describes the spacetime behind the shock-wave 
is consistent with the fact that the galaxies appear to be uni- 
formly expanding. In this paper we construct such a model 
assuming critical expansion, (k = 0), for the FRW metric, and 
what emerges is a new, essentially exact solution of the Einstein 
equations. We show that reasonable physical requirements on 
the TOV equation of state put an interesting constraint on the 
possible position of the shock-wave relative to the center of the 
explosion. Using this constraint, we derive precise estimates for 
the shock position at present time, as predicted by this model; 
that is, at the time in this model at which the Hubble constant 
and the value of the background radiation temperature agree 
with observed values. The constraint on the shock position can 
be interpreted as a new length scale that is derived from the 
model, and this length scale is not determined by any adjustable 
parameters in the problem other than the experimentally deter- 
mined values of the Hubble constant and the background radi- 
ation temperature. The constraint on the position of the inter- 
face arises because the interface is a true shock-wave. There 
is no similar constraint on the position of the interface in the 
well-known Oppenheimer-Snyder model, where the interface is 
a contact discontinuity, |H| [T3f . 

In this paper we start with a critically expanding, (k = 0), 
FRW metric under the assumption that the equation of state 
agrees with the equation of state that applies in the standard 
model of cosmology after the time of the thermal uncoupling of 
matter from radiation. This uncoupling occured at a tempera- 
ture of approximately 4000°i^, at about 300, 000 years after the 
Big Bang in the standard model, flU, [T9|1 . In Section we derive 
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a system of ODE's that determine the TOV metrics that match 
the given FRW metric across a shock-wave interface, (equations 
( 15731) , Q5751) below). In fact, we derive the shock equations in the 
case of a general FRW metric, allowing for k ^ and for general 
equations of state, (equations ( |3.41| ), ( |3.42| ) below). In Section 
H we give a derivation of the FRW equation of state in terms of 
the cosmological scale factor. This includes a discussion of the 
FRW metric in the presence of both matter and radiation fields, 
assuming that the pressure due to matter is negligible, and that 
there is no thermal coupling between the fields. The results in 
Section |] also apply for arbitrary k. 

To obtain the ODE's for the TOV metric, we must rework the 
theory in [[13], [T5[] where a given outer TOV metric is the starting 
point, instead of a given inner FRW metric which we require 
here. These ODE's, which are non-autonomous, simultaneously 
describe the TOV pressure p and the FRW shock position r, 
assuming conservation of energy and momentum and no delta 
function sources at the shock. (We let barred quantities refer to 
TOV variables and unbarred quantities to FRW variables, c.f. 
|}I74j.) We then derive a formula for the TOV energy density p, 
(the only remaining undetermined variable in the TOV metric), 
which, together with solutions of the ODE's, determine the TOV 
solutions that match the given FRW metric across a shock-wave 
interface. The ODE's take a particularly simple form when the 
cosmological scale factor R of the FRW metric is taken to be 
the independent variable instead of the TOV radial variable f, 
the independent variable in the usual formulation of the TOV 
system, |T8|, . In Section [5| we present a rather complete phase 
plane analysis of these equations and we prove that there exists 
a unique bounded orbit. This orbit describes the TOV pressure, 
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but does not constrain either the initial shock position or the 
TOV energy density. We show that along this orbit, the pressure 
jump across the shock-wave has the property that the ratio of 
the TOV pressure to the FRW pressure at the shock is equal to 
a slowly varying function that is bounded between 1/9 ~ .1111 
and a = y/17 — 4 « .1231, where the FRW pressure is supplied 
by the background radiation. Using this bound on the TOV 
pressure along the orbit, we obtain the following sharp upper 
and lower bounds for the (squared) distance that the shock- 
wave can propogate over and above the (geodesic) motion of 
the galaxies, as a function of "starting time", R*. (Here R = 1 
denotes present time in the model, and we view the starting time 
R* < 1 as the earliest time at which the shock-wave solution 
has settled down to the point where our model applies; that is, 
as entropy increases, we expect shock-wave solutions to settle 
down to simple time asymptotic configurations, and we assume 
here that this time asymptotic solution agrees with our model 
from R* onward.) The inequality reads, (see |7.31|) below), 

(2.62 x <r ,_ r2< (2.65 x KTQTg lp ( ±_ \ 



KqHq \R^j HqHq \R* 

Here the distance r is given in terms of the Hubble length 

.98 

Hq 1 « — x 10 10 lightyears, 
h 

where 

Hq = 100/io km sec' 1 mpc' 1 , 
and it is generally agreed that ho lies in the interval .5 < ho < 
.85, M. For example, if we take T = 2.736 w 2.7°K, R* = 



6 



2.7/4000, and ho = .55, (a recently quoted value), the above 
estimate reduces to 



2 2 

r — r\ 



/.019\ 2 
\Ho~) ' 

We conclude that the distance the shock-wave has traveled, 
(over and above the motion of the galaxies), between R = = 
2.7/4000 and present time R = 1, as predicted by this model, is 
approximately .019 times the Hubble length. 

In the standard interpretation of the FRW metric in Cosmol- 
ogy, the galaxies are in freefall, and traverse geodesies r = const. 
Thus we can interpret r 2 — r 2 above as the (squared) distance 
that the shock-wave travels over and above the motion due to 
freefall, a result of the fact that mass and momentum are driven 
across the shock-wave as it evolves outward. From this point of 
view it is a bit surprising that the quantity r 2 — r 2 is independent 
of the starting position r*. 

Using the formula for the TOV energy density, we next prove 
that the minimal physical requirement p > p > p > 0, then 
places an additional contraint on the initial shock position r* 
that depends on the starting time R*. We prove that once this 
constraint is met at one time, it is met at all succeeding times 
in the solution, and the density and pressure profiles are phys- 
ically reasonable. Putting these results together, we obtain the 
following upper and lower bounds on the shock position at the 
present time as a function of background radiation temperature 
Tq, the Hubble constant Hq, and the value of the scale factor R* 
at which we start the shock-wave, (c.f. ( |/.5Y|) and (|Y.38|) below): 
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r > ff -{(5.1xl0-)gf(| 



r < m 



\ 



1 + 



.76 

TO 



fl. + (2.6X10-^ in (i-) 



The maximum shock position is plotted in Figure 2 for the case 
Tq = 2.7°K, and ho = .55. In Section § we compare these bounds 
to the analagous bounds one obtains in the case of pure radia- 
tion, thus making contact with the exact solution discussed in 
P|. For example, at Tq = 2.7° K we obtain 



36/i 36/W1 + 2.5i? t 
< r < 



Hq ~ ~ Hq 
c.f. ( |ST7|) below. (Of course, since we are neglecting the matter 
field, we do not have R/R = Ho at the same time when T = To 
in the pure radiation model.) We note that in this case the 
distance from the shock position to the center of the explosion 
is significantly beyond the Hubble length. 

In summary, starting with the idea that there might be a 
shock-wave that marks the outer boundary of the expansion 
that we measure by the Hubble constant, one's first reaction is 
that nothing quantitative could be said about the position of 
the shock without knowing details concerning the nature of the 
spacetime beyond the shock-wave, or details about the mecha- 
nism that might have created such an explosion in the first place. 
And to a large extent this must be true. But what we find inter- 
esting here is that this simplest shock-wave cosmological model, 
consistent with both the observed values of the Hubble constant 
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and the background radiation temperature, contains within it 
unexpected constraints on the possible position of such a shock- 
wave, and the shock position is comparable to the Hubble length. 
That is, we find it intriguing that this length scale comes out 
of the model, and provides an answer that apriori needn't have 
been so reasonable. 

In conclusion, we ask whether our expanding universe could 
have evolved from the center of a great explosion that gener- 
ated a shock-wave at its leading edge. If so it makes sense to 
wonder whether some of the far away objects that we observe 
in the nightime sky are possibly due to similar explosions that 
originated at other locations in spacetime. We now know that 
the scale of supernovae is not the largest scale on which classical 
explosions have occured in the universe. Indeed, it was reported 
in a recent issue of Nature, that on May 7, 1998, a gamma ray 
explosion emanating from a faint galaxy known as GRB971214 
erupted, and for two seconds the burst was more luminous than 
the rest of the universe combined. This is the largest explosion 
ever recorded, and redshifts place it at about 12 billion lightyears 
away. Moreover, conditions at the explosion were equivalent to 
those one millisecond after the Big Bang in the standard model. 
Thus we pose the question: could explosions such as this, or 
even greater than this, have given rise to our own "expanding 
universe"? Indeed, could we then observe other similar explo- 
sions in distant regions of spacetime beyond the expansion of 
our own universe, (that is, beyond the shock-wave that marks 
the edge of the expansion we measure by the Hubble constant)? 
We propose the shock-wave model presented in this paper as a 
natural and simple starting point for a futher investigation of 
these issues. But independently of this, the model provides a 
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new, essentially exact solution of the Einstein equations that we 
feel is interesting in its own right. 



2 Preliminaries 

According to Einstein's theory of general relativity, H, the grav- 
itational field is described by a Lorenzian metric g that satisfies 
the Einstein equations 

G = 8 -fr, (2.1) 

on 4-dimensional spacetime. Here G is the Einstein curvature 
tensor, Q denotes Newton's gravitational constant, c denotes the 
speed of light, and T is the stress energy tensor, the source of the 
gravitational field. In this paper we are concerned with FRW 
and TOV metrics, two spherically symmetric metrics which are 
exact solutions of (|27l| ) when T takes the form of a stress tensor 



for a perfect fluid, namely 

Tij = (p + pc 2 )u t Uj + pgij, (2.2) 

where p denotes the mass-energy density, p the pressure and 
i, j = 0, 3 denote indices of spacetime coordinates. The FRW 
metric is given by 

ds 2 = -diet) 2 + R 2 (t) I l -^dr 2 + r 2 dU 2 ) , (2.3) 

{ 1 — kr l J 

and the TOV metric is given by 

ds 2 = -B{f)d{ct) 2 + A(r)- l df 2 + f 2 dQ 2 , (2.4) 
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where dfl 2 denotes the standard metric on the unit two-sphere. 
We write the TOV metric in barred coordinates so that it can be 
distinguished from the unbarred FRW coordinates when we do 
the matching of these two metrics below, c.f. |T3j. Assuming co- 
moving coordinates and substituting ()2.3j) into the (|2.1|) yields 



the following FRW equations, [g O, |TB|, [TO 



R 2 = ^pR 2 ~ k, (2.5) 



and 



R P so a\ 

P=-P~ M - (2-6) 

The unknowns R, p and p in the FRW equations are assumed 
to functions of the FRW time t alone, and "dot" denotes dif- 
ferentiation with respect to t. Assuming co-moving coordinates 
and substituting (|2.4|) into the ( |2TT|) yields the following TOV 
equations, 

— = 47rf 2 A (2.7) 
dr 

and 

B' , s 

— = -2-^— , 2.9 

where 

A = l-™. (2.10) 

c z r 
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Here the unknown functions are the density p, the pressure p, 
and the total mass M, [jT2|, which are assumed to be functions 



of f alone, and prime denotes differentiation with respect to f. 
In the next section we fix an FRW metric and derive equations 
for the TOV metrics that match the given FRW metric across a 
shock-wave interface at which the metric is only Lipschitz con- 
tinuous, and across which conservation of mass and momentum 
hold, and at which there are no delta function sources. In |T3j 
it is shown that the shock surface is given implicitly by 

4-7T 

M(f) = -p{t)f\ (2.11) 



and the metrics (|2T3Q , flZTgj ) are identified via a coordinate trans- 
formation in which 

f=Rr. (2.12) 



3 Derivation of Equations 

In this section we derive equations that describe the time evo- 
lution of an outgoing spherical shock-wave interface together 
with an outer TOV metric, such that the shock surface matches 
a given FRW metric on the inside, and such that conservation 
of energy and momentum hold across the interface. The main 
point here is that we are assuming a given inner FRW metric, 
rather than assuming a given outer TOV metric as in [|T3| , |T5 



Thus we seek a pair of equations that determine an outer TOV 
metric that matches a given FRW across a shock-wave interface. 
Rather than deriving the shock equations, we shall write them 
down and prove that solutions of these equations determine a 
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shock-wave solution of the Einstein equations. (The reader can 
obtain a formal derivation of these equations by reversing the 
steps in the arguments below.) 

Equation ( pj.lUj) in our first theorem below is the first equation 
in the pair of ODE's that we will work with. 

Theorem 1 Assume that p(t), p(t), and R(t) solve the FRW 
system 



R 



8tt£ 
\ 3c 4 



R 2 p - k, 



(3.1) 



over some interval 



Assume that 



and that 



P = -3^(p + p), 



I = (t h t 2 ). 



R(t) > 0, 



(3.2) 



(3.3) 



(3.4) 



i?^0, (3.5) 

on I. We assume WLOG, (by the choice of positive square root 
m 07H)), that 



R>0. (3.6) 

Assume further that r(t) is a positive invertible function defined 
on /, and define f(t) on I by 
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r = Rr. 

Define functions M(f) and p(r)) by 



47T 

M(f(t)) = —p(t)f(t) 



(3.7) 



(3.8) 



and 



p(r) = 



M'(f) 

4:7Tf 2 



(3.9) 



where prime denotes differentiation with respect to r. Assume, 
finally, that r(t) satisfies 

1 f p-p \ l-kr 2 

R \p + p) Rr ' [ ' } 

for some function p, and that p, p, M, p, and p are all positive 
valued functions on I. Then for all t € I we have, 



p 



where, /p~dj], 



= 



A = 1 



lOp-p 
1-70 ' 

A 

1 — &T 2 ' 

2QM 



(3.11) 



c 4 r 



and 



7 = 



P + P 
P + P' 



(3.12) 
(3.13) 

(3.14) 
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That is, Theorem 1 implies that for a given FRW solution, ( |3.1U|) 
implies the conservation condition ( |3.11|) when M and p are 
defined by (|3~8| ) , and (|3T^), (these latter two equations being the 
shock surface matching condition and the second TOV equation, 
respectively, [[13], |T5|). Here dot denotes ^, and we assume c = 1. 

Proof: Differentiating ( |3~5|) with respect to ct and using 
gives 



But ( ETB1) gives 



»V dM - * — 2- 

M = ——r = 4:TTpr r. 
dr 

47T 



M 



/6r + 47rpr r, 



so from (p.lh| ) and (|XTT3|) we get 



Rr 



Hp- p) 



p- 



Using (|3T2| ) in (\SA{\ ) gives 



(3.15) 



(3.16) 



(3.17) 



y P + P 

r = —rti— . 

P~ P 

Using (|3.Y| ) and simplifying we have 



(3.18) 



fR + rR 



-Rr 



P- P> 



Using ( |3.1(J|) to eliminate f from Q3.19|) gives 



(3.19) 



1 — kr 
R 2 r 2 

We now use the identity 



'P + P\ (P + P^ 
,p - p) \p — V> 



(3.20) 
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1 — kr z 



1 



R 2 r 2 1-0' 
which follows from (|3.1j) and (|3~8| ) . Indeed, 



• 9 2Q M „ 9 , 
i? 2 = -?r^rR 2 - k. 



c 4 r 3 



But ( g7T3| ) implies that 



and using this gives 



-^r- = C 1 " A )?> 



1 - A , 



or, (c.f. 10 0), 



(3.21) 



r 2 i? 2 = -A + fl -kr 



Using in (gggp gives (frZIp , as claimed. 

Now using ( |3.21| ) in (|3.2U|) yields 



1 



'p + /A f p + p" 



(3.22) 



1 — \p — pJ \p — p J 
Solving this for p gives ( |3.11|) , where we have used ( |3.14j ). This 
completes the proof of Theorem 1. 

For a given FRW metric, Theorem 1 tells us that the ODE 
( |3.1U| ) can be taken in place of the conservation constraint ( |3.11|) , 
and the reversal of the steps in the above proof can be regarded 
as a formal derivation of the ODE ( |3.1U|) . For later convenience, 
we now record the following additional equations that follow 
from the hypotheses of Theorem p.l| . 
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Corollary 1 Assume that the hypotheses fa. 1\ ) through fa. 1 Q) 
of Theorem \3. 1\ hold. Then the following equations are valid: 



p + p 7# 



p — p 7$ — 1 ' 
6 



r 



(p + p) = Vl - kr 2 (-^==^j (p-p). 



e 



'p + v\ fp + P^ 



1 — 6 \p — p) \p — p) ' 

f = — ^— Vl - kr 2 ^/T^6, 
7# — 1 

r 2 R 2 = -A + (1 - kr 2 ) , 



1 — kr' 



1 



r 



1-0' 



(3.23) 

(3.24) 

(3.25) 

(3.26) 
(3.27) 

(3.28) 



Proof: By Theorem 4.1, we know that ( |3.11| ) holds, and using 
this in the LHS of ( gggj) gives the RHS of ( gggj ). Also, from 
(^22j), 



ri? = Vl - kr 2 Vl-0, 
and using this in ( |3.19|) gives 



VP - p/ 



(3.29) 



(3.30) 



Using (ftZ3p in flCTp gives (frZBp . From (|XTTj) we get 
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p — p = 



lOp-p p(l-j0) (p + p)9) - {p + p) 



so 



1-70 1-^6 



J9 — p 6 — 1 



1-70 



/? + J9 1-70 

To verify ( [j.24|) , we use ( |3.2(j| ) which we write in the form 



(3.31) 



r 



Vl — kr 2 



1 



1-6 



VT^lO-l 



6. 



and so from ( |3.31|) we have 



r 



Vl — kr 2 fp — p^ 



(3.32) 



y/1 — 6 \p — p, 
Solving for (p + p)f in ( [j.32|) gives ( ft.24|) . Finally, to obtain 
( ggg) , equate the RHS's of flXgUfl and flX52p . Equations ( gg7p 
and (|3.28|) have already been derived as (|3.22|) and ( |3.21|) within 
the proof of Theorem |T|. This completes the proof of the Corol- 
lary. 

Now assume that p(t), p(t), and R(t) solve the FRW system 
( |371|) and ( |3T2| ) for t <E I, and assume that the hypotheses ( pT3|) 
to (p.l(J| ) of Theorem [T] hold. We know from Theorem [I] that 
the conservation condition ( ft.ll|) also holds. We now find an 
equation for p(f), (equation ( fj.at>|) below), which guarantees that 
p solves the TOV equation (pTS|) , since then, in light of (|3T9|), the 
functions p(f), p(f), and M(f) will then solve the TOV system 
as well. Defining A(f) by ( |3.13| ), we can define the function B{f) 
as a solution of the ODE, [T3], |T5|, 
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B'(r) 



2p'(f) 



B(r) p{r) + p(r) ' 
thus determining a TOV metric of the form 



(3.33) 



ds 2 = -B{f)dt 2 + A-\f)df 2 + f 2 dQ 2 , (3.34) 
that solves the Einstein equation G = kT for a perfect fluid with 



stress tensor 



Tij =V9ij + (p + p)UiU 3 . 

For this metric, co-moving coordinates are assumed, fflPfl , and 
thus the 4-velocity u is given by 



w u = V5, u % = 0, i= 1,2,3. 

Note that we are free to choose any positive initial value for B 
by suitable rescaling of the time coordinate i. The next lemma 
demonstrates that if p satisfies equation ( |3.56|) below, then as a 
consequence it also satisfies 



V 



QMr - ( \ 47ror 3 \ , , 



c 4 r 2 



M 



(3.35) 



which is equivalent to the TOV equation (|2.8|) 



Lemma 1 The hypotheses ( \3. 1\ ) to 1 u\ ) of Theorem \3. 1\ , to- 
gether with the equation 



p = 



QM (I - kr< 



A 



c 4 r 2 



Rr 



1 — kr 2 , 



(p-p) [l + ^A-\ (3.36) 



imply that p(r) also solves the TOV equation 
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Proof: By we have 



VI — kr 2 



1 



Rr VT^e 7 

and substituting this together with (|3.12Q into ( |3.3(j| ) gives 



V 



QM 

c^f 2 



\/l — kr 2 



But using ( |3.24|) and (|3.8|) in (|3.3Y|) we obtain 



P = 4-2 



^ f( , +p - ) f 1 + ^! u - 1 



(3.38) 



which directly implies the TOV equation (^TB|). This completes 
the proof of the Lemma 0} 

Our results now imply the following theorem which intro- 
duces the system of ODE's whose solutions we analyze in sub- 
sequent sections: 

Theorem 2 Assume that p(t), p(t), and R(t) satisfy the FRW 
equations 1\) and fort G /, and that the other hypotheses 
(fTTSp through (i^.lUj) of Theorem \3. 1\ hold. Assume further that 
(r(t),p(t)) solves the system of ODE's 



p 



1 (p — p\ 1 — kr" 



R Kp + pJ Rr 

gm(p-p) (1 + 3J 



c 4 r 2 



Rr 



(3.39) 



(3.40) 
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for t G /, where r, M(r), and p(r) are defined for r(t\) < r < 
f(t 2 ) by $![]$) and Then p(f), p(f), M(f), solve the 

TO V system '/[ ), ([%. and the conservation condition ^3. 1 1\ ) 
holds for all t E I. Furthermore, under these assumptions, the 
system is equivalent to the system 



dr 
dR 

dP_ 
dR 



P-P\ (I- kr r 



R(QR 2 - k) \Q + P f 

_1 (Q + 3PKP-P) 
2 QR 2 - k 



r 



for R(t{) < R < R(t 2 ), where 

(Q,P,P) = ^(p,p,P), 

has the dimensions of inverse length squared. 



(3.41) 
(3.42) 



(3.43) 



Note that the equivalence of system ( |3.39| ) , ( |3.4U|) with ( |3.41|) , 
( |3.42| ) follows because of the assumption R ^ 0. This also implies 
that P and Q can be considered as functions of i?, in which 
case equations (|3.41|) and (|3.42|) closes to form a well-defined 
nonlinear system of two ODE's in the unknowns r and P. After 
solving ( |3.41|) , ( p.42|) the dependence of R on t can be recovered 
from (|3.1|) . Thus for a given FRW metric and a given solution 
of (|3.41|) , ( |3.42|) , the only variable remaining to be determined 
is the TOV energy density Q on the outside of the (outgoing) 
shock-wave. To obtain a closed form expression for Q, write the 
shock surface equation ( |3.8| ) in the form 



2QM = Qf c 



(3.44) 
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But the second TOV equation for M is 



dr 

which we can rewrite in the form 



dM 9 . 

= 4irpf 2 , (3.45) 



^-(2QM) = SQf 2 . (3.46) 
dr 



Differentiating ( |3.44j ), substituting into ( |3.4t)|) and solving for Q 
yields the formula 

fl-is^-^TS"- (3 ' 47) 

Note that if Q decreases as the shock moves outward, (that is, 
the f position of the shock increases), then the second term 
in ( p7T7|) is negative, and so Q < Q, (the density behind the 
shock is greater than the density in front of the shock), as is 
the case for classical shock-waves in fluids, [T7J. Note, however 
that the physically necessary condition Q > 0, or the physically 
reasonable condition Q > P, is not guaranteed, and depends on 
the particular solution. 

The final theorem of this section tells us that solutions of the 
ODE's ( |3.4(J| ), (|3.41|) do indeed determine exact shock-wave so- 



lutions of the Einstein equations when a (suitable) FRW metric 
is given. 

Theorem 3 Assume that p(t), p(t), and R(t) satisfy the FRW 
equations 1\ ) and ( \3. "4) for t E /, and that the hypotheses 
C f^. through of Theorem [O] hold. Assume further that 



(r(R), P(R)) solve the system of ODE's ( \FJT\) , for Rjh) < 



R < R(t2). Assume that Q, P, M, Q, P, and A are all positive 



22 



and that the shock speed in FRW coordinates is less than the 
speed of light throughout the interval I. Then there exists a C ' 
invertible coordinate transformation mapping (t, r) — > (t, f) of 
the form 



t = t{t,r), (3.48) 
f = f(t,r) = R{t)r, (3.49) 

such that, under this identification, the resulting TOV metric 
matches the given FRW metric Lipschitz continuously across the 
shock surface r = r(t). (The angular coordinates 6 and (j) are 
implicitly identified.) Moreover, the Lipschitz continuous metric 
defined by taking the FRW metric for r < r(t) and the TOV 
metric for r > r(t) defines a shock-wave solution of the Einstein 
equations fl\Q, \Qj), c.f. j73j [T5[/ . In particular, the Rankine- 
Hugoniot jump conditions 

[T tj y = 0, i = 0,...,3, (3.50) 

hold across the shock surface; there are no "delta function sources '' 
on the surface; there exists a regular C > coordinate transfor- 
mation defined in a neighborhood of each point on the shock 
such that the metric components in the transformed coordinates, 
(which can be taken to be Gaussian normal coordinates), have 
smoothness level C ' 1 ; and the matched metric determines a weak 
solution of the Einstein equations in the sense of the theory of 
distributions, c.f. fiTSjJ. 

Proof: The existence of the coordinate transformation is proved 
T3j pages 278-280 under the assumption that the shock sur- 



m 



face is nowhere characteristic in the sense of (4.43) of that paper. 
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Using (4.55) of the same reference, the non-characteristic con- 
dition can be re-written as 



r 7^ 



A 



(3.51) 



which holds here because we assume that ^4 > 0, r > 0, R > 
and f > 0. Since the normal vector n to the shock surface is 
non-null, (because we assume that the shock speed is less than 
the speed of light), and the functions c(t, r) and c(t, f) in Lemma 
9 of [|T3f are here equal to R(t)r and f, respectively, it follows 
that the conclusions of Lemma 9, JT3J are valid. Moreover, the 
conservation condition ( |3.11|) is valid, and thus the argument 
in |T3J that leads to (3.9) in this latter reference, implies that, 
under our hypotheses, condition (5.5) of Lemma 9, follows 
from the conservation condition ( |3.11| ) above. (Note that the 
condition (2.20) of |T3j], assumed in that paper, is not needed 
here.) Since the conclusions of Theorem |3| are just a re-statement 
of the conclusions of Lemma 9, |jI3fl , the proof of Theorem |3| is 
now complete. 

A remarkable aspect of the formulation of the shock equations 
given in (|3.41|) and ( |3.42| ) is that, if Q and P are given functions 
of R, (which can be obtained from the FRW equations once an 
equation of state is specified), then the equation ( |3.42|) for P 
uncouples from the f equation ( |3.41|) . Thus, in principle, one 
can solve system ( |3.41| ) and ( |3.42|) by first solving the scalar 
non-autonomous equation (|3.42|) for P, and then plugging the 
solution P(R) into ( p.41|) to obtain a scalar non- autonomous 
ODE for the shock position r. 

As an application, and to clarify the way system ( fj.41|) and 
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(|3.42|) works, we now demonstrate that the set of shock-wave 
solutions of the Einstein equations determined by system ( p. 41 



(|3.42|) includes, as a special case, the exact solutions first pre- 
sented in 0. To verify this, set k = in flOTl ) and flcT42D to 
obtain the system 



dr 
dR 

dP 
dR 



1 



P-P 




QR 3 \Q + P t 

1(Q + 3P)(P-P) 
2 Qi? 



(3.52) 
(3.53) 



Now it is easy to verify that the solution in [|P|, (as given in 
equations (5.1) through (5.14) of that paper), satisfies the fol- 
lowing equations: 

P = crQ, (3.54) 



for some constants < a, a < 1, and 

-3(l+cr) 



l+3cr 



(3.55) 



Q 


f- 


Qo 






-( 







(3.56) 



(3.57) 



where the subscript zero denotes values at some particular time, 
say present time in the FRW solution. For our purposes here, 
we now take ( |5.54|) to ( |3.hY| ) as an ansatz, and show that this 
ansatz is consistent with system ( |3.b2|) , ( |3.b3|) , and that system 
( |3.b2|) and ( |3.b3| ) then determines the other relations in (5.1) to 
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(5.14) of |I|J, (including the relation between a and a). To this 
end, set Rq = lf\ and c = 1, and substitute (|5.b4| )- (|5.bt)| ) into 
equation (|5.b^| ) to obtain 



RdQ 3(l + a){a-a) 
QdR ~ ~2 ~d 



(3.58) 



But from ( p.btiQ , 



RdQ dln(Q) 



= -3(1 + a), 



(3.59) 



QdR d\n(R) 

and putting this into the left hand side of ( |3.b8|) and solving for 
a gives 



a 



a(7 + a) 



3 1 



(3.60) 



which is the relation between a and a given in (|4.^|) -(|4.^|) of \\T4 
Substituting ( |i5.54j ) and ( |3.55|) into equation and sim- 

plifying gives 



1 dr A 



1 



'3cr — a s 



2dR QR 2 \3 + aJ' ^' 6 ^ 
and using ( |3.b6| ) and ( |3.bY| ) in ( |3.61|) gives, after simplification, 



3d 



u 



1 + 3<t) (3 + ^) 



(3.62) 



3 Note that we are free to fix the initial condition Rq = 1 because the FRW metric 
is invariant under rescalings of R when k = 0. Indeed, under the scaling R — > aR, the 
variables Q, P, P and f are in varia nt, so th e shock position r = f/R rescales like r — > a r. 
It is easy to see that system ( 3.52 ), ( 3.53 ) is invariant under this rescaling. 
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Equation ( |3.b^ ) determines the constant tq from the constants 
Qq and a or a. To connect this with the formulas in fTfljl , note 
that (|5.bt>|) and ( |5.bY| ) together with the matching condition 



r = Rr, 



yeild the formula 



Q = Qo - 

\r , 



(3.63) 



(3.64) 



Since we take i?o = 1, 



Q 



Qo4 



(3.65) 



But ( |3.65|) checks with jI4j because (3.2)-(3.4) there imply that 

(3.66) 



_ 8ttQ 87r^7 

Q = —z-p = 



r- 



where, (using the notation of |T4|), 

1 / 

7 = 



(7 



2tt£ V1 + 60- + CT 2 , 
A calculation using ( |3T0D| ) verifies that 



(3.67) 



3(7 — 0" 



(3.68) 



(l + 3a)(3 + a) 

and thus ( |3.62| ) implies (|3.66|) as claimed. 

We can now solve for the density Q using the TOV equation 

dM 



dr 



(3.69) 



together with the shock matching condition 
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M = — -pf 3 . 
3 P 



That is, using ( |3.6h|) in ( |3TT0| ) gives 



M(r) = 47T7r, 



and thus (|3.t)U|) implies 



P 



1 



which implies that 



and 



P = <jQ, 



(3.70) 



(3.71) 



(3.72) 



(3.73) 



(3.74) 



also agreeing with the formulas arrived at in [ptj] . Finally, sub- 
stituting ( |3T56|) into the FRW equation 



and integrating gives the time dependence of as 



(3.75) 



[2 + 3a)VQ~o 



2+3<r 



(3.76) 



where we have used the initial condition R(0) = 0. 

For our construction below, the case a = | is relevant, (the 
case of pure radiation). In this case, 



a = V17-4w .1231..., 



(3.77) 
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as follows from the formula 



a = l-V9a 2 + 54a + 49 - |a - -. (3.78) 

Since this is the only case relevant to the discussion below, we 
will from here on let a denote the special value a = y/YJ — 4. 



4 The Equation of State 

In this section we derive the FRW equation of state that we use 
to close the system ( |3.35|) and ( |3.3(j| ). In this paper we consider 



the case when the equation of state for the FRW metric agrees 
with the equation of state in the standard model of cosmology 
after the time that the radiation in the universe uncoupled from 
thermal equilibrium with matter. This is an effort to account for 
the observed microwave backgound radiation level in our shock- 
wave model. Our idea is that if the expanding universe arose 
from a great explosion, then one might conjecture that the ex- 
pansion would have settled down to a uniform expansion by the 
time that this decoupling occured. In the standard model of cos- 
mology, the thermal uncoupling of radiation and matter occured 
at about 300, 000 years after the Big Bang, at a temperature of 
about 4000 degrees Kelvin, []]], [TTj] . Thus, we analyze our 
shock-wave model in the case that there is an energy density 
p r for the radiation, (which supplies a pressure p r = (l/3)p r 
via the Stefan-Boltzmann law), and a separate energy density 
p m for the matter, which is assumed to exert a zero pressure 
p m = 0. Since the scale factor R is the independent variable in 
our shock equations ( |3.35|) and (|3.3b|) , we now obtain formulas 



for the FRW energy density and pressure as functions of R. 
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We start with the FRW equations in the form, |I9j| , 



and 



d 

dR 



(4.1) 



(4.2) 



(Again, we assume that "dot" denotes d/d(ct).) We can rewrite 
O as 



R 2 = QR 2 - k, 



and equation ( |4.2|) as 



_d_ 
dR 



(QR 3 ) = -3PR 2 



(4.3) 



(4.4) 



where Q = ^fp and P = ^fp have dimensions of inverse length 
squared, c.f. (pT2T|) . Now assume that the energy in the FRW 
system is in the form of pure radiation and matter alone, so that 



Q = Qr + Qm (4.5) 

where Q r , Q m denote the (appropriatly scaled) energy density 
of radiation and matter, respectively. Further, assume that the 
pressure of radiation is given by the Stefan-Boltzmann Law, []T9| 



Pr = (l/3)Q r , 



(4.6) 



and that 



P m = 0, 



(4.7) 
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so that the uncoupling implies that 



P = Pr + Pm = Pr = (l/3)Q r . (4.8) 

Finally, assume that 

Qm = (4.9) 

for some positive constant (3, so that the total energy of matter 
within a (geodesically) expanding volume in the FRW metric 
remains constant. Substituting (f|15j) through into (|T4j) 

gives 



j(Qr# 3 ) = 1 

which has the solution 



(4.10) 



Qr = ^ (4-11) 

for some positive constant a. Then the FRW pressure is given 
by 

P = Pr = (4-12) 

We conclude that the equation of state that applies to the FRW 
system under the assumption that radiation is uncoupled from 
matter is given by 

P = Jf- (4.14) 
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Putting ( |4.l3|) , ( |4. 14|) into equations (|3.41|) , (|3.42|) gives the sys- 
tem of ODE's 



dr R (a- PR 4 \ (1 - kr 2 ) 

dR ~ (3a + 0R- kR 2 ) \3a + (3R + PR 4 ) r ' 

dP _ 1 (3a + PR + 3PR 4 )(a — PR 4 ) 

dR ~ ~2R5 3a + PR - kR 2 ' ^ ' ' 

Observe, again, that a nice feature of the formulation ( |4. 15j) , 



( f4. 16| ) is that, in this formulation, the second equation ( |4.16| ) for 



the TOV pressure P uncouples from the first equation ( |4.15| ) for 
the shock position r. 

Equations ( |4.13| ) and Q4.14Q together with the Stefan-Boltzmann 
law imply that the temperature of radiation is proportional to 
1/R. Indeed, let T = T(R) denote the temperature of radiation. 
The Stefan-Bolzmann law relates the energy density of radiation 
p r to the temperature through the relation 



p r = aT 4 , (4.17) 



where 



a « 7.664 x 1Q~ 15 ^-(K°)\ (4.18) 

cm 6 

Since Q r = ^rfp r , we can write this as 



Q r = aT 4 , 



where, IT9 



a = 



SirQa 
3^ 



(4.19) 



(4.20) 
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defines the constant a. Now by (|4.11| ) 



3a 



= Qr 



8tt£ 



Pi 



™ 3c 4 " 
we have the following lemma: 



8tt£ 

3c 4 



aT 4 = aT\ 



Lemma 2 The Stefan- Boltzmann law implies that 

W 1/4 I 

~ ItJ r' 



(4.21) 



5 Restriction to k = — Phase Plane Analysis 



We now analyze system ( |4.15| ), (|4. in the case of critical ex- 
pansion when k = 0. The case k ^ will be considered in a 
subsequent paper. To start, recall that in the case k = 0, sys- 
tem ( |3.41| ), (|3.42|) reduces to the (non-autonomous) system 



where again 



dr 
dR 

dP_ 
dR 



'P-P* 



r 



(QRZ) \Q + P / 

1(Q + 3P)(P-P) 
2 QP 



(Q,P,Q,P) = -g-j-tp.p.p.p), 



(5.1) 
(5.2) 



(5.3) 



all have the dimensions of inverse length squared. Assuming 
now that Q and P are given by ( |4.l3|) , Q4.14|) , 



Q = — + ^ 
W R A P 3 ' 
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p = ^ 



(which models a universe of matter and radiation assuming no 
thermal coupling), and substituting this into ( p.'Zj) yields 

the system 



dr R ( a- PR A \ 1 



dR (?>a + (3R) \3a + /3R + PR 4 J r 
dP 1 {3a + (3R + 3PR 4 ){a- PR 4 ) 



(5.4) 



(5.5) 



dR 2R 5 3a + (3R 

which is just system (|4.15|) , (|4. in the case k = 0. Solutions 



of system (|5~T4] ) and ( |5T5| ) determine the shock position r(R), 
(the position as measured by the radial coordinate of the FRW 
metric that is behind the shock-wave), together with the TOV 
pressure P(R) in front of the shock. The TOV pressure p(f) is 
then recovered from the solution (r(R), P(R)) by inverting the 
equation f = Rr(R) and using P = ^rp- The function R(t) is 
obtained by solving the FRW equation (|X1Q with k = 0, and 



due to the scaling law for this equation, we are free to choose 
the scale factor Rq such that Rq = 1 at present time in the 
universe. The constants a and (3 that determine Q rescale with 
choice of Rq, and are determined from initial conditions for the 
FRW metric. For an FRW metric that models the expanding 
universe, we can take one of the two initial conditions as 

Q = 3a + f3 = Hi (5.6) 

where H is the present value of the Hubble constant, c.f. ( [[TTJ ) . 
Finally, the TOV energy density Q is given by the formula in 
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Q 



:(Qr 3 ) = Q + ——Q. 



(5.7) 



3r 2 dr ~" ' 3 dr 

This simplifies under our special assumption (|4.l5|) for Q. In- 
deed, by ( FCT3Q , 



df di? G?f 



d (3a \ ( dr" 
dR\R* + l$)[ r+ dR, 



Since 



dr d(Rr) _ 
— - = r + R — 



where "' 



(5.8) 



(5.9) 



dR dR dR' 
dR is given by equation (|5TT|). Putting (|5TT|) and ( p~8|) 



into ( |5.Y|) and simplifying yields the following expression for the 
TOV energy density Q : 



Q = Q 
where 



(4a + 0R)(3 + ^R + w){3 + ^R)ar 2 
R A {ar 2 {3 + ^R + w)(3 + Jr) + (1 - w)R 2 } 



(5.10) 



PR 4 



w = 



a 



(5.11) 



We conclude that each choice of constants a and (3 and each 
choice of initial conditions for ( |5.4j ) and ( |5.5|) determines a shock- 
wave solution of the Einstein equations, at each point where all 
the variables are positive. We now analyze solutions of system 
(gTp and O in detail. 

Substituting w for P in (|b.4Q , (|b.b|) yields the equivalent sys- 
tem 
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dr 
dR 

dw 
dR 



R 

1 



w 



(3 + gi2) \3 + £R + w 



ar 



4w ( (1-w)(3 + £R+3iu)' 
~R\ 8(3 + gii> 



(5.12) 



(5.13) 



Now ( |5.13| ) is a non- autonomous scalar equation for P that un- 
couples from equation f |5.12|) . In order to do a phase plane anal- 
ysis of equation (|5.13|) , and to analyze the behavior of solutions 
as R — > oo, we now rewrite (|5.13| ) as an autonomous system. To 
this end, set 



S = — Rj 

a 



and 



u 



1 

S' 



(5.14) 



(5.15) 



Substituting these into ( |5.13| ) and letting "dot" denote d/dS, we 
obtain the following autonomous system of two ODE's that is 
equivalent to equation ( |5.13| ): 



du 2 



dw 

W — — — = 4:WU 

dS 



1 



1 - w)(/3 + 3(1 +w)u) 
8(/? + 3w)w 



(5.16) 



(5.17) 



We now analyze the phase plane associated with system ( |5.1t)|) 
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System ( p.lY| ) has a line of rest points at u = 0, and 

an isocline where w = 0. Setting the RHS of equation ( |b.lY|) 
equal to zero gives 



1 



[1 -w)(fl + 3(1 + w)u) 



8(f3 + 3u)w 
and solving this for w gives 



0. 



(5.18) 



w 



(3/3 + 8u) 
2u 



■1 + 



1 , 4 (f3 + 3u) 

-L- — ii 

\ 3(3/? + 8u) 2 



(f>(u). (5.19) 



Thus the isocline is defined for < u < oo by 



Note first that when (3 
isocline degenerates to 



w = (j)(u). (5.20) 
0, (the case of pure radiation), the 

<j)(u) = Vl7-4 = a. (5.21) 

It is straightforward to verify that when (3 = 0, the isocline is 
also a solution orbit of system ( |5.16|) , ( |5.17| ), and the special 
solution in [jT3j| , discussed above starting with (|3.54|) - (|3.5Y|) , cor- 
responds to this orbit. The special value a = Vl7 — 4 also is 
important in the case (3 ^ 0. The next theorem gives the quali- 
tative behavior of the solution orbits of system ( |5.1(j| ), (|5.1Y| ) in 
the (w, w)-plane when (3 ^ 0. 

Theorem 4 Assume that (3^0. Then the following statements 
are true regarding the phase plane of system (jp-lOj) , (jo. 
(Here we define an orbit of system ( \b. lb\ ), ( \o.l'Jj ) to be a function 
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w = w(u) such that (u(S), w(l/S)), (u = 1/S), is a solution of 
system ( \h~TQ ), ($7T1\)): 

(i) The isocline w = <j)(u) defined by (|5.iyQ is monotone increas- 
ing for < u < oo, and satisfies 



Jim0(w) = u = v 17 — 4 ~ .1231..., (5.22) 
lirn^O) = i « .1111..., (5.23) 
lim^W = ~i ^ -01097... ( 4) 



fuj Orbits can only cross the isocline w = <f){u) once, from right 
to left in the (u,w)-plane, as S increases, (see Fig 1.). 

(Hi) Along any orbit w = w(u) we have 

lim w(u) = o. (5.25) 



(iv) There exists a unique orbit w cr i t (u) satisfying 

lim w crit (u) = ~. (5.26) 

Moreover, all orbits w = w(u) starting from initial conditions 
(uo,wq) such that wq > 4>(uq), (that is, starting above the iso- 
cline), satisfy 

lim w(u) = oo; (5.27) 

u— >0 

and all orbits starting from initial conditions (uq, wo) such that 
wo < 4>(uo), (starting below the isocline), satisfy 

lim w(u) = — oo. (5.28) 
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Proof: To verify (|h>.22|) , we have 



lim 6{u) = 



r (3/? + 8u) 
lim 

r (3/3 + 8w) 
lim 



.1 + . i + +!").„ 



\ 3 (3/3 + 8w)< 



•1 + 



\ 



1 + 



16 



= -1 + \Zl7 = (7. 



We next show that </>(u) tends to a monotonically from below 
as u — > oo. Note that by (|5.19| ), w = <j)(u) is equivalent to 



0= 1 

which we rewrite as 



(1 -w)(P + 3(1 + w)u) 
8(/3 + 3w)w 



(5.29) 



= 8(/3 + Su)w -(l-w)[(3 + 3(1 + w)u] . (5.30) 
Now differentiating (|5.3U|) implicitly with respect to u gives 



8(/3 + 3w) — + 24w = - W + 3(1 + w)u] — 
du du 



3 + 3w + 3u 



dw 
du 



(5.31) 



Simplifying we obtain 



(9/3 + 24u + 6uw) 



dw 
du 



-3(w 2 + 8w - 1). 



(5.32) 



Now the roots of w 2 + 8w — 1 are 
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ct = V17-4w.1231, o = -\/l7-4, (5.33) 
and thus we conclude that, along the isocline w; = 0(it), 



^ > if iy < a, (5.34) 
aw 

— < if w > a, 5.35) 
du 

where we use the fact that 4>{u) > for all u. Thus, it suffices to 
show that (f)(u) ^ o for any u in order to conclude that 4>'(u) ^ 
for < u < oo. So assume for contradiction that (j)(u) = a. But 
solving for u in ( |b.5U|) gives 



(itf - cr)(w + |(j|) V ' 

and thus w = a leads to a contradiction unless (3 = 0. We 
conclude that if /5 7^ 0, then 0(w) monotonically increases to a 
as w — > 00, thus proving (|5.22| ). Statement ( |5.^3| ) follows from 
( |5.3U| ), and (|5.23|) follows from (|5.34 ). Thus the proof of (i) is 



complete. 

Statement (ii) follows because w = only on w = 0(it), 
w; > if w > 0(w), and since we have shown that 4>'(u) > 0, it 
follows that orbits can only cross the isocline from right to left 
in forward (S-time. 

To verify (iii), we show that all orbits tend in backward time, 
(increasing u), to w = a. To see this note that 



lim w = lim 4wu 



8(3 + » 
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4wu < 1 



w 



8w 



(5.37) 



where approximately means to leading order as u — * oo. Now 
each orbit that starts above w = (j){u) decreases as u increases 
unless the orbit crosses the isocline, in which case the orbit in- 
creases from there on out as u — > oo. It follows that orbits start- 
ing below w = 4>(u) can never cross w = <fi(u) at any value of u 
larger than the initial value. Thus, since lim^oo 4>{u) = <r, all 
orbits must be bounded above in w by the maximum of {a, wq} , 
and bounded below by the minimum of {1/9, Wq} . But from 
( |b.3Y|) , we must have that 



lim 

u— >oo 



1 



1 



W 



8w 



= 0. 



(5.38) 



Indeed, if not, then (|5.3Y| ) implies that \w\ tends to infinity as 
u — >■ oo, which implies that w is not bounded as u — > oo, and 
this contradicts the above bounds. Since o is the only positive 
root of \ 1 — \ , we conclude from ( |5.38| ) that 



lim w = a. 

u— >oo 



This completes the proof of (iii). 

We now give the proof of (iv). Our approach for this is to 
write the ODE for w as a function of u along orbits, and study 
the limit u — > 0. From (|5.1b|) and (|5.1Y|) , 



dw 
du 



4w 
u 

4w 
u 



1- w)(/3 + 3(1+ w)u) 



S(P + 3u)w 



w 



8w 



(5.39) 
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where approximate equality means to leading order as u —> 0. 
Now assume for contradiction that there exists an orbit w = 
f(u) that is bounded in a neighborhood of u = 0, but such 
that lmv->o f(u) ^ 1/9. The boundedness condition implies that 
( |5.39| ) applies with errors that are bounded as u — > 0. That is, 



~ = -(l-if^+OW (5-40) 
du u \ 8w J w v ' 

9w — 1 

it 

where 0(1) denotes a constant that depends on the bounds for 
w but is independent of u as u — > 0. Integrating ( |5.41| ) leads to 
the estimate 



9w ~ 1 = f?f + (l)^W, (5.42) 
9u>o — 1 V w / 



where (uq, wq) are taken as initial data, uq > 0. But (|b.42| ) im- 



plies that if lim u _,o^ 1/9, then w = f(u) is unbounded near 
u = 0. From this we conclude that every orbit that is bounded 
as u — > satisfies 



lim it; = 

We now show that there exists at least one orbit such that 
lim M ^o w = 1/9. Note first that any orbit starting from initial 
data (uo,Wq) that lies on the isocline, wq = (p(uo), uq > 0, 
must lie above the isocline for all < u < uq because we know 
that 4^ < on this interval, and (f/(u) > 0. Since the isocline 



(5.43) 
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decreases to o « .1231 > 1/9 as u — > 0, it follows that wq > 1/9 
for initial data lying above the isocline, and hence lim u ^o w > 
1/9 along an orbit starting from such initial data. But our 
argument above shows that when this happens, we must have 
lim^o^ = +oo. We conclude that lim^o^ = +oc for any 
orbit starting from initial data above the isocline, wo > 4>(uo). 
Similarly, if the initial data (uq, wq) lies below the line w = 1/9, 
that is, wo < 1/9, then also wq < (p(uo) because we have that 
(j){u) > 1/9. Thus from ( ggg ), g > 0, and so it follows that 
lim u _^o < 1/9, and our argument above implies that lim u ^o w = 
— oo. We conclude that lim^o^ = — oo for any orbit starting 
from initial data below the line w = 1/9; and lim^o w = +oo 
for any orbit starting from initial data above the isocline, wq > 
4>{uq). Now consider all orbits emanating from initial data on 
some fixed vertical line u = e > 0. Then if wq > 0(e), we have 
lim u ^o = +oo; and if w$ < 1/9, we have lim u _>o = — oo. So 
define 



where the limit is taken along the orbit emanating from the 
point (e, wo). We now claim that the critical orbit emanating 
from initial condition (e, wo) satisfies lim^o^ = 1/9. To see 
this note first that w + > 1/9 because orbits below w = 1/9 tend 
to —oo as u — > 0. We show next that the orbit emanating from 
(e, w + ) cannot tend to w = +oo as u — > 0. To see this, note that 
if lim^o^ = +°° along the crititcal orbit, then this must be 
true for all orbits starting in a neighborhood of (e, wq) as well. 
Indeed, if lim^o^ = then at some positive value of u we 
must have w > <f>{u) along the critical orbit; and so by conti- 




(5.44) 
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nuity, nearby orbits must also rise above the isocline at some 
u > 0, and hence by above we know that lim u „»o w = +oc along 
orbits sufficiently close to the critical orbit. But this contradicts 
the fact that W- is a greatest lower bound. We conclude that we 
cannot have lim M ^o w = +°° along the critical orbit. Similarly, 
we cannot have UtOu-^qW = — 00 along the critical orbit because 
then nearby orbits would also satisfy lim w ^o w = —00 since they 
would cross w = 1/9 before u = 0, and again this would con- 
tradict the fact that w + is a greatest lower bound. Since we 
cannot have lim u ^o w - -00 or lim w ^o w = +00, it follows from 
( |5.4;j| ) that the only alternative is that lim M _>o if = 1/9 along the 
critical orbit, as claimed. 

We now show that the critical orbit is unique. To this end, 
rewrite equation ( |5.1Y|) as 



w = 4wu 



[l-w){/3 + 3{l + w)u) 



8(/3 + 3w)w; 

Differentiating ( |5.45|) with respect to w gives 



= F(u, w) 



dw dF 



dw dw 



= 4u 



1 + 



6w; + j3/u 



> 0. 



(5.45) 



(5.46) 



8(3 + P/u) 

But ( |5.4ti| ) implies that the distance between orbits is increasing 
in forward time S, (that is, increasing as u = 1/5 decreases. 
Indeed, 



dF 

[w 2 - W\) = F(u, w 2 ) - F(u, Wi) = — (li, W*)(W2 - Wi) > 

ow 

(5.47) 
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if W2 — W1 > 0. This implies that there cannot be two orbits that 
satisfy lim^o^ = 1/9 since the distance between them would 
then tend to zero as u — > 0, contradicting (|b.4V|) . This finishes 
the proof of (iv), and thus the proof of the theorem is complete. 



The salient properties of the phase plane for system Qb.lb|) , 
( |5.1Y| ) are sketched in Figure 1. Note that as (3 — > 0, the iso- 
cline moves up to the line w = a, (continuously, except for a 
jump from 1/9 to a at u = 0, (3 = 0). The isocline is a curve 
of absolute minima of orbits that cross the isocline, and the 
isocline, together with all orbits, tend to w = a as u — > 00, 
R — > 00. Moreover, all orbits except the critical orbit tend to 
infinity as u — > 0, (R—> 00), and so the critical orbit is the only 
orbit bounded for all values of R > 0. Along both the critical 
orbit and the isocline, the following apriori bounds hold for all 
< R < 00 : 



1/9 w .1111 < w < u « .1231. (5.48) 

Note, however, that the critical orbit and the isocline do not 
coincide except in the limiting case (3 = 0, in which case both 
reduce to the line w = a, which also can be indentified with 
the special solution constructed in In particular, Figure 

1 describes how this special solution is imbedded in the larger 
class of solutions that allow for general initial data. 



6 Conditions for Q > and Q > P 

In this section we obtain conditions which guarantee that Q > 
and Q > P, physically reasonable conditions on the TOV energy 
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density and pressure, (Q = ^rp~ where p is true TOV energy 
density.) In particular, Q > P guarantees that Q > whenever 
the solution orbits of system ( |5.16|) , ( |5.1Y| ) satisfy w = PR 4 /a > 
0. We begin with the formula (|5.1U| ) for Q : 

Q-Q + P R )$ + gg + ^)( 3 + foV' (61) 

W V m{ar 2 (3 + ^R + w)(3 + ^R) + (l-w)R 2 }' K ' J 



Asking that the RHS of ( |5.1(J|) be positive, and using the formula 



Q = — + ^ 
we see that Q > is equivalent to 

a \ar 2 (3 + -R + u/)(3 + -fl) + (1 - w)i? 2 ) (6.2) 
[a a ) 

3 

-Ua + (3R) (3 + -R + w)ar 2 > 0. 
a 

Solving flET21) for r 2 leads to the following inequality that is equiv- 
alent to Q > : 



ar 2 < - — a—^ . 6.3 

3 + 2-R + w 



Equation (|6.3|) implies that the condition Q > puts a contraint 
on the maximum possible shock position at a given value of R. 
The following theorem implies that if the condition holds at 
some value R = R* in a solution of ( |5.16|) , (|5.1Y|) , then it holds 



for all R > R* in that solution, so long as < w < 1 and 
dw/dR < 0. Both of these conditions are satisfied along the 
critical orbit where 1/9 < w < a « .1231. 
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Lemma 3 Define the quantity {}/ by 

(I - w)R 2 



3 + £R + iw 



ar 



Then for any solution of ( \b.l v 4 ), ( \q.1 j| j we have 



(6.4) 



at each point where 



and 



d_ 
dR 



{}/>0 



< w < cr, 



(6.5) 



(6.6) 



dw 
dR 



< 0. 



(6.7) 



Lemma |3| implies that {}/ is monotone increasing along any 
solution of ( |5.16| ), ( |5.1Y| ) that satisfies (|677|), and thus if ( |B~3|) 
holds at a value R = It** in such a solution, then it must hold at 
all R > FL. 



Proof: Starting with ( |6.4| ) we have 



_d_ 
dR 



2(1 -w)R 



?>+ [ -R + w 



a 



dr 2 n d 
+ R 2 



w 



2(1 -w)R 



dR 
2(1 -w)R 



dR [3 + ILr + Wjij 



(6.8) 



3 + ^R + w (3 + ^ + ^(3 + ^) 



II 



where we have used ( p.l4 ). This simplifies to 



47 



_d_ 

dR 

Moreover, 



2(1 -w)R (2 + £R\ +R2 _d^ {} 



2, + ^-R + w \3 + £R 

a \ a 



di2 



77- 



(6.9) 



{}// = 



1 — it; 



di? 1 3 + 2-R + iy 
{Z + 2-R + w) 



7/ 



(!-»)(! + 



/3 _|_ riu> ■ 



(3 + ^ + ^) 2 



_ (3 + gfi+l)(-j)-(l-«>) 



3 + &R + iy 



> 



(i-«05 



(3 + ^.R + u; 

V a 



> 0. 



(6.10) 



where we have used ( |6TB|) and (|677| ). Using ( |6.1U| ) in (|6.9|) and 
simplifying gives 



d {h > ; (1 : w) ( 4 +f^ 5 > o. 



3 + ffl + io)) (3 + ffl 



(6.11) 



This completes the proof of Lemma |3|. 

We now obtain a corresponding condition for Q > P. Using 
( |h.ll| ) and ( |5TB|) we know 



- aa 



(6.12) 



Using this together with the formula ( |b.lU|) for Q, it follows that 
the condition Q > P will hold if 
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R 4 Q 



3 + fR 



a 



v a ' 3+^R+w 



R?(l-w) ii J - °- 



Solving (|b.l5|) for ar gives the equivalent condition 



1 + cr 



-ar 2 < {} //7 , 



where 



{} 



in 



R\l- w) 
3 + + 



Thus to get ( |b.l4| ) it suffices to have 



^«r 2 =(l + e)ar 2 <{} /7/ , 



where 



4a 



3-a 

We conclude that Q > P holds so long as 

{{}/// " (1 + e)«r 2 } > 0. 
Lemma 4 Define the quantity 6?/ 

{}iv = {{}m-(l + e)ar 2 }. 
Then for any solution of ( \h.l°4) , (\o.l$) we have 



d 

dR 



{}iv > 
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at each point where ( \6.t\ ) and fib. '/[ ) hold. 



In particular, Lemma [| implies that if (|b.l8|) holds at a point 
R* in a solution of (|h.l2|) , such that ( |BT5|) and ( |6TT| ) hold 

for all R> R*, then we conclude that (|OU]) holds at all points 
R > R*, and thus that Q > P for all R> R*. 



Proof: Differentiating we obtain 



d 2(1 — w)R , . dr 2 « d 
X}iv = - Br, ( 1 + e )«^ + i? 



w 



dR lj % + Pr + w v ' dR dR[3 + £R + w. n 

a v. ' a ' 11 

2(1 -w)R 2(1 -w)R 2 d 

- 3 + ^ + W _(1 + e, (3 + ^ + »)(3 + ^) +fl d7j {} " 

2a-^/2- £ + f ^ +fll i 



Now using ( p.lU|) in Qb.21|) and simplfying yields 



±{}iv > ^[ g 3 + ^fl(4 + gfl - 2.) > 0. (6.22) 
ait 3 + £ i?+w a a 

a 

since 2e < 4 and k; < 1. This concludes the proof of Lemma [|. 
We have proven the following theorem: 

Theorem 5 Assume that fib. b\) and fib. '/[) hold for all R > R* 
on a solution of fib.1%) , fib.lty) . Then {}/ > at R = R* is 



equivalent to Q > at R = it**, and implies that Q > for all 
R > R*; and if {}iv > at R = R*, then we must have Q > P 
for all R > R*. The condition {}/ > is equivalent to 
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2 (l-w)R 2 

ar 2 < -± — ^ , 6.23 

3 + ^R + w V ; 



and the condition {}iv > simplifies to 



2 /l-cr/3\ (l-w)R 2 . . 

ar 2 < - '— v R ' . 6.24 



7 Estimates for the Shock Position 



In this section we take system ( |5.12|) , (|5.13| ) as a simple cos- 



mological model in which the FRW metric behind the shock- 
wave at position r is assumed to model the expanding universe. 
Given this, we now estimate the position of the shock-wave in 
the present universe as determined by this model. In this model, 
the expanding universe is modeled by an FRW, (k = 0,) met- 
ric in which the energy density Q and pressure P are given by 
( ^TT3| ), (^23), that is, the same as that assumed in the standard 



cosmological model after the time of thermal decoupling of mat- 
ter with radiation, (approximately 300,000 years after the Big 
Bang in the standard model, [01|). The FRW metric is assumed 



to have been created behind a radially expanding shock-wave 
due to a great explosion into a static, spherically symmetric 
universe modeled by a TOV metric. Given these assumptions, 
we have shown that conservation of energy at the shock then 
implies that the position r of the shock-wave is determined by 
equation ( |5.12|) , where r is the radial coordinate in the FRW 



universe behind the shock. Equation (|5.12| ) is coupled to equa- 
tion Qb.lai) for the TOV pressure P, and the TOV energy density 



Q is then given by the formula ( |h.l(J|) . In this section we assume 
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that w = R 4 P/a lies on the critical orbit w = w cr i t (S), (S = ^). 
(This is justified by the fact that, according to Theorem f|, this 
is the only orbit bounded for all R, and all orbits are asymp- 
totic to this one as R — > 0.) By Theorem |], w ranges between 
a and 1/9 along the critical orbit, and thus we have the apriori 
estimate 



1/9 « .111K w < a w .1231. (7.1) 

The only remaining piece of information missing is the initial 
condition for the shock-wave. At first one might think that this 
initial condition can be chosen arbitrarily, but as we have shown 
in the last section, the condition that the energy density be 
positive in front of the shock-wave, or that it be larger than the 
pressure in front of the shock, puts a constraint on the maximum 
shock position at a given time. That is, assuming that w lies 
on the critical orbit implies that the hypotheses of Theorem |5] 
hold, and thus condition ( |6.23|) , 

2 (l-w)R 2 
3 + -R + w 

a 



is equivalent to Q > 0, and the condition ( |b.24|) , 



2 l-a3\ (l-w)R 2 
ar 2 < '— ^ — ^ , 7.3) 

is sufficient to guarantee that Q > P, at any given value of 
R. Moreover, if ( |b.23| ) or ( |b.24j ) hold at a given value R = R*, 



Theorem [5] tells us that they continue to hold for all R > R*. 

Under the above assumptions, we now obtain estimates for 
the shock position. To start, rewrite ( |5.12|) as 
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dr A 



2(1 -w)R 



dR a(3 + lR)(3 + £R + w) 
2(1 -w)fR 



R + %(3 + w)) [R + 3% 
Using (|7TT| ) in (fTTSj ) gives the estimate 



(7.4) 



dr 2 
< < 



(l-wJ)%R 



(fl+(3 + iu + )f) (R + 3%) ' dR (fl+(3 + w_)||)(i2 + 3g 



(7.5) 



where w_ = 1/9 « .1111 < w + = a « .1231. That is, 



2(l-a)i2 



dr 2 a 
< < 



(16/9)fl 



P (R + (3 + *)§) (fl + 3§) ^ ^ 2 (fl + (f? ) (R + 3§ 



Now by direct calculation, the solution to the ODE 



P 



(7.6) 



dr z 



= C- 



R 



dR {R + A){R + B)' 
for positive constants A, £>, and C, is given by 



(7.7) 



r = + In 



^ , AC „ ^ -BC 

R + A\^B ( R + B\^B 



,R* + AJ \R* + B, 



(7.8) 
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where inequalities can be substituted for equalities in ( [777]) , 
(|T%1) . Applying this to (|775Q gives the inequalities 



7% > In 



1 R + A. 
M* + A 



+ ■ 



R + B + 



[7.9) 



where 



r 



r„ < In 



' R + A 
M* + A^ 



■ R + B 

M^TbZ, 



(7.10) 



A 
B 



a = 



b = 



2(3 + w)(l 



6(1 — w) a 



w f3 2 ' 



(7.11) 



and t4_, ^4 + , are obtained by substituting W-, w + for w, respec- 
tively, in the above expressions, etc. 

We now evaluate a and (3 in terms of the present value of the 
Hubble constant Hq and the observed microwave background ra- 
diation temperature Tq. Here we let subscript zero denote value 
at present time in the FRW metric, and WLOG we assume that 
Rq = 1. Recall that the FRW equation ( |2T5| ) for k = can be 
written as 
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R 



= QR' 



so that 



#o — vQo, 

where Qo denotes the present value of the (scaled) energy density 
in the universe at present time. By ( |4.13|) , 

Qo = 3a + (3, 

where 3a is the energy density of radiation at present time, and 
(3 is the energy density of matter at present time. Let T = T(R) 
denote the temperature of radiation. Then ( fQTp is 



(7.12) 



where (|4.2U|) gives 



a = 



3^' 



Setting Rq = 1 and solving ( |4.2l|) for a gives 



a, 
3 



a = o T o 4 5 



(7.13) 



and using this in ( |4.13|) gives 



= Q Q - 3a = H£ - aT 4 . (7.14) 
We evaluate the above constants using the values, ([[1911), 
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f = 7.425 x 1CT 29 era g~\ (7.15) 

c = 2.997925 x 10 10 cm sec' 1 , (7.16) 

Ity = 9.4605 x 10 17 era, (7.17) 

mpc = 10 6 pc = 3.2615 x 10 6 Ity. (7.18) 

a = 7.5641 x 1(T 15 erg cm' 3 K \ (7.19) 

H = 100/io km sec' 1 mpc' 1 , (7.20) 

T = 2.736 °K, (7.21) 



Here, Q is Newton's gravitational constant, c the speed of light, 
Ity is lightyear, rape is megaparcec, °K is degrees Kelvin, a is the 
Stefan-Boltzmann constant, Tq is the observed microwave back- 
ground radiation temperature |TTj, and Hq is Hubble's constant, 
where ho is generally accepted to be between .5 and unity. (We 
take ho w .55 as a recently quoted value.) Using these values we 
calculate 

a = 4.6852 x 10" 27 Ity' 2 K'\ (7.22) 

H = 1.023/io x 10~ 10 Ity' 1 . (7.23) 
Using the above values we obtain from ( |Y.13| ) and (|7TT4J) that 
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= (1.492 X 10- 7 )^| I = 8 - 34 X 1 2 " 6 , (7-26) 
where we used << 1 at the approximate equality. Using 

■"0 

these values we can evaluate: 



- = (4.66 -^-^ r ° 
P 



A + = (3 + .1231)- = (4.66 x KT 7 )-|, (7.27) 







2 ' 



5 + = 3- = (4.48 x 10-^ 

a + - 2(3 + .1231)(1 - .1231) 
a+ ~ .1231 ^ - (b.wy x iu 

6(1 - .1231) a , fix T n 4 

6+ = — ^ — — = (6.377 x 10 _(h 



.1231 (3 2 v ' 'hlHl 1 

A = (3 + l/9)^ = (4.64xl0- 7 )§, (7.28) 
p ho 



B- = 3- = (4.48 x 10-0-|, 
a = ^±1Z^M « = (7 . 427 x 1(r * ^ 



1/9 v ' y ^ 2 ' 



6(1 -1/9) a ^ 1co ^ 1n _ 6 , T ( 



6_ = 777^-=; = (7.162 X 10 



1/9 p v * 7 ^ 2 - 
Now, assuming that the uncoupling of matter and radiation 
occured at a temperature less than 4000 degrees Kelvin, [|I9| , it 
follows from ( |4.^1| ) that 
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R* > 2.2/4000 = 6.75 x 10~ 4 , (7.29) 

and so it follows that we can essentially neglect the A's and B's 
in estimates ( |7T^|) and fl7TTDp , given their small values in ( |/.2Y|) 
and ( |/.28|) , and assuming this, estimates (|7!3| ) and ( |<M0|) reduce 
to, 



(o+ + 6+) ln(l/R,) <r 2 - rl < (a_ + 6_) ln(l/i?*). (7.30) 

Using ( |Y.2Y|) and ( |Y.28|) to estimate ( |Y.5U| ) gives the following 
estimate for the distance the shock-wave must have traveled 
between r = and R = 1 as predicted by our model: 



; 2 .62 x io- W ln m < r2 _ < (2.65 x : io-^ ln / 1 > 



JIqHq \R*J ^Hq \R 

(7.31) 

Here the distance r is given in terms of the Hubble length 

#o xlO 10 . (7.32) 

In particular, ( |7.31| ) shows that, in this shock-wave model, the 
quantity r 2 —r\ is essentially independent of the starting position 

7V 

As an example, if we take ho = .55, Tq = 2.736°K and R* = 
2.7/4000 in ( |Y.31| ), we obtain the estimate 



2 2 

r — r ^ 



'0 

In the standard interpretation of the FRW metric in Cosmology, 
the galaxies are in freefall, and traverse geodesies r = Const. 
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Thus we can interpret r 2 — r\ in ( |Y.3l| ) as the (squared) dis- 
tance that the shock-wave travels over and above the motion 
due to freefall, a result of the fact that mass and momentum 
are driven across the shock-wave as it evolves outward. We 
conclude that the distance the shock-wave has traveled, (over 
and above freefall), between R = = 2.7/4000 and R = 1, 
as predicted by this model, is approximately .019 of the Hubble 
length. (Recall that f = R(t)r measures distance in lightyears 
for the three dimensional space at fixed time t in the k = FRW 
metric.) 

We now discuss the initial condition r = at R = R*. We 
saw in ( |tj.24| ) that the condition Q > P put constraints on the 
maximal shock position at each value of R. Using the value 
a = .1231 in (|b.24j ) gives the inequality 



.759^ 



.759 



R* 



(3.11)a + PR i + 



(4.64xlQ- 7 )r 4 



m 



(7.34) 



Estimate ( |7.34j ) is the bound on the initial shock position, im- 
posed by Q > P, in terms of the Hubble length. Putting ( [7.34Q 
together with ( |7.31| ), we conclude that the maximal distance 
r max from the shock-wave to the center of the explosion r = 
at present time R = 1, given as a function of starting time i?*, 
2.7/4000 < R*l, (assuming the shock-wave started at position 
r = r* at R* > 2.7/4000, and such that is restricted by ( |^.34|) 
so that Q > P for all R > R*), is predicted by this model to be 
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r, 



H, 



-l 



\ 



1 + 



.76 

(4.6xlQ- 7 )T4" 



i^ + (2.6x 10- 7 )^|ln (-?-) 

tin \K*J 



(7.35) 

For example, taking the value ho = .55 and Tq = 2.736°K gives 
the formula 



\ 



.76 



1 + 



(8.5x10-5) 

m 



i?* + (4.9x 10- 5 ) In (^-)- (7-36) 



This function is plotted in Figure 2. Putting (|Y.35|) together 
with ( |/.31|) we obtain the following upper and lower bounds for 
the shock position r at present time R = 1 assuming that it 
starts at R = R*, and such that Q > P holds for all R> R* : 



J- n 



r > H^l^.l x HT 4 )^- 



In ( 



h\ \R 



(7.37) 



r < H Q ~ l 



\ 



1 + 



76 T 4 / 1 \ 



(7.38) 



8 The Case of Pure Radiation, (3 = 

As a point of comparison, in this section we redo the calculation 
of the shock position under the assumption (3 = in Q4.l3|) ; 
that is, under the assumption that the energy density Q is due 
entirely to radiation. (See and the solution discussed at the 
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end Section |3].) Thus assume that a = aT 4 /3 is as given in 
(|7TT3|), but that (3 = 0. We estimate the position of the shock- 
wave in this model at the time R = 1, where T = Tq. Now of 
course, since (3 is determined in ( |7.14j ) from Hq in the above 
analysis, the value of ^ in the pure radiation model will not 
coincide with Hq at the time when T = Tq. Nevertheless, for 
comparison purposes, we shall estimate the radial position of 
the shock-wave in the pure radiation model at time R = 1 in 
terms of the Hubble length Hq 1 given in ( |/.32|) . 



In the case (3 = 0, the constraint ( |b.24|) that guarantees Q > 
P reduces to 

\ (l + ff)(3 + iu) / v ' 

and the critical orbit becomes w = a. Using w = a ~ .1231 in 
|8TTC|) gives 

49 

r, < ^i?,. (8.2) 



(Note that in the alternative case a = 0, the case of pure matter, 
the RHS of tends to infinity, and thus Qtj.24| ) places no 



constraint on the shock position. This is consistent with the 
fact that when a = 0, the pressure is zero, and the shock-wave 
reduces to a contact discontinuity. For example, Q = 0, P = 0, 
solves the shock equations ( |3.41|) , (|5.42|) and it is not difficult 
to show that the solution of the shock equations in this case 
reduces to the k = version of the Oppenheimer-Snyder model, 
first presented in [|T4|. In these Oppenheimer-Snyder models, 



there are no constraints on the shock position corresponding to 
(£23).) 
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Setting (3 = and w = a in (|b.l2|) gives 

dr 2 2(1 - a) R 



(8.3) 



dR 3(3 + a) a 

as the differential equation for the shock position. Integrating 
gives 



r 2 = 



1 — a) R 2 2 

f h r*. 

3(3 + a) a 



(8.4) 



Using ( |8.2| ) for the maximum value of yields the following 
bounds on the shock position r at the time R = 1 when T = Tq 
that are analagous to (|V.3V|) and ( |7T3S| ) and apply when (3 = : 



1 



1-cr) 1 

< r < 



1-^ 



yfa\ 3(3 + o)~ ~ ^\ 3(3 + a) 

I From (|7TT3|) and (5^3 ) it follows that 

1 

118/io^o" 1 - 



+ .24iL 



(8.5) 



a 



(8.6) 



Using this value together with the value a = .1231 in (|875| ) yields 

36/i . . 36/ioVl + 2.5R* (Q ^ 

-JT < r < 77 • ( 8 - 7 ) 

Note that the shock position at R = 1 that applies to the exact 
solution given in [jnj], which was discussed in detail at the end 
of Section is the case R* = in flST7j) . 
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Figure 1: 
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